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摘 要 : 参照 《中 国 植被 》 中 的 植被 分 类 体系 ,结合 野外 考察 结果 ,建立 了 适合 中 国 西北 农 牧 交错 带 


的 植被 分 类 体系 。 以 覆盖 研究 区 的 多 幅 Landsat 影 
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现 了 对 研究 区 植被 信息 的 提取 。 提 取 时 , 先 利用 完全 约束 的 最 小 二 乘 模型 对 遥感 影像 进行 混合 像 
元 分 解 ,将 整个 研究 区 划分 为 植被 区 和 非 植被 区 ;在 植被 区 ,基于 光谱 特征 、 纹 理 特征 和 地 形 特征 ， 
构建 CART 决策 树 ,获得 了 乔木 林 、 灌 从 和 草原 等 7 种 主要 植被 型 组 ;在 植被 型 组 内 ,基于 不 同 植被 
类 型 NDVI 的 季节 差异 特征 ,构建 NDVI 差 值 比值 指数 (NDVI_DR) ,将 乔木 林 和 灌 从 区 分 为 常 绿 和 和 
落叶 植被 型 ,使 用 温度 植被 干旱 指数 (TVD1) ,将 草原 划分 为 荒漠 草原 、 典 型 草原 和 草 甸 草原 3 种 类 


型 ,从 而 得 到 各 个 植被 型 的 空间 分 布 范 围 。 经 验 i 


下 ,最 终 分 类 的 总 体 精 度 能 达到 79.51% ,kappa 系 


数 为 0.773。 本 文 所 采用 的 分 类 方法 充分 利用 了 遥感 数据 既 有 的 光谱 信息 和 纹理 信息 ,同时 辅 以 


地 形 信 息 。 实 践 结 


表明 ,分 层 分 类 和 多 种 指标 相 结合 的 方法 可 以 有 效 实现 对 影像 跨 幅 的 、 以 复杂 


镶嵌 结构 为 主要 特征 的 农 牧 交错 带 植被 信息 提取 ,精度 较 高 ,技术 可 行 。 
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植被 的 研究 一 直 以 来 都 备 受 关注 "| ,因为 掌握 
植被 的 分 布 特征 有 利于 推动 生态 环境 建设 ,协调 区 
域 可 持续 性 发 展 。 本 文 所 选 的 研究 区 西北 农 牧 交错 
带 , 常 年 降水 较 少 、 莹 发 量 大 且 植 被 稀 朴 ,是 典型 的 
生态 脆弱 区 。 

在 西北 农 牧 交错 带 , 由 于 地 形 和 地 貌 的 复杂 性 、 
气候 的 特殊 性 以 及 人 类 活动 的 影响 ,植被 类 型 的 变 
化 极为 复杂 。 以 往 的 研究 表明 ,对 于 西北 农 牧 交错 
带 这 样 的 地 表 覆 盖 复 杂 区 域 ,采用 “分 步骤 ,分 层 
次 ,分 类 别 ” 的 方法 能 有 效 的 提高 最 终 的 分 类 精 
度 2 ] 。 此 外 , 主 成 分 变换 、 织 帽 变换 (又 称 K-T 变 
He) 波段 优选 和 一 些 常 见 的 遥感 指数 的 组 合 也 可 
以 提高 分 类 精度 "2 Khatami R 等 ”发现 纹理 
言 息 和 多 时 相 图 像 的 使 用 能 够 改善 分 类 的 整体 精 
度 。 遥 感 影像 能 够 捕获 植被 生长 变化 中 的 一 些 生 理 
特征 和 物候 规律 ,而 这 些 特征 大 多 具有 季节 性 ， 
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因此 使 用 这 些 季节 性 变化 的 特征 并 结合 多 时 相 的 遥 
感 影像 能 够 充分 利用 遥感 影像 反映 出 的 特征 完成 植 
Ba REI! 。 在 具体 操作 时 ,一般 采 用 不 同时 相 
NDVI 的 组 合 来 提取 植被 的 落叶 和 常 绿 信息 "1]。 上 
述 遥 感 植 被 分 类 方法 ,多 集中 于 植被 覆盖 度 较 高 的 
区 域 ,但 是 将 以 上 方法 结合 起 来 用 于 西北 农 牧 交错 
带 植 被 信息 提取 的 研究 较 少 。 

本 文 以 遥感 影像 为 主要 数据 源 ,探索 出 适合 于 
地 表 履 盖 复 杂 植被 稀疏 区 域 的 植被 信息 提取 方法 。 
本 研究 最 终 得 到 的 西北 农 牧 交 错 带 的 植被 分 布 为 深 
入 分 析 该 区 域 地 表 水 热 格局 提供 基础 数据 ,为 该 区 
的 生态 保护 提供 理论 文 撑 。 


1 研究 区 概况 与 方法 


1.1 研究 区 概况 
农 牧 交错 带 是 一 种 植被 类 型 .气候 特征 和 生产 
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方式 均 具 有 过 渡 性 的 特殊 地 理 区 域 。 本 文 所 研究 的 
西北 农 牧 交错 带 地 理 位 置 为 105°35' -110°54'E， 
36°49’ -40°11'N, 位 于 鄂尔多斯 高 原 和 毛乌素 沙 地 
区 域 ( 如 图 1) ,为 干旱 区 和 半 干 旱地 区 。 西 北 农 牧 
交错 带 具 有 典型 的 生态 脆弱 性 ,其 南部 为 毛乌素 沙 
地 、 北 临 库 布 其 沙漠 、 西 接 乌 兰 布 和 沙漠 。 研 究 区 海 
拔 范 围 在 640 ~2 060 m 之 间 , 总 体 呈 现 出 西高 东 
低 、 南 高 北 低 的 特点 。 研 究 区 总 面积 约 达 8 x 10° 
km ,包含 内 蒙古 .宁夏 陕西 的 10 个 县 ( 旗 ) 。 植 被 
类 型 主要 为 适合 干旱 和 半 和 干旱 环境 的 乔木 灌木 和 
人 工 栽 培植 物 。 
1.2 分 类 体系 

本 文 以 《中 国 植被 ) 中 植被 分 类 体系 "为 参 
考 ,充分 利用 可 获取 的 研究 区 其 他 相关 资料 和 多 次 
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实地 考察 所 获取 的 数据 制定 了 相对 合适 
系 。 其 中 考虑 了 以 下 原则 : 

(1) 以 优势 种 为 主 :植被 类 型 特征 的 分 类 标准 
由 植物 群落 中 的 优势 种 确定 。 对 于 本 研究 区 ,植被 
类 型 多 为 荒漠 植被 ,不 可 能 包含 所 有 《中 国 植被 》 中 
提 到 的 植被 类 型 。 因 此 将 研究 区 内 的 植被 亚 型 直接 
分 为 植被 型 ,将 研究 区 分 布 极 少 或 无 分 布 的 类 型 不 
列 人 本 分 类 体系 。 

(2) 考虑 动态 特征 :考虑 植被 的 次 生性 质 和 演 
替 关 系 , 将 部 分 植被 型 合并 。 例 如 ,不 再 将 灌 草 丛 单 
独 划 分 为 一 个 植被 型 ,而 是 将 其 纳入 灌 从 所 属 的 植 
被 型 中 。 

(3) 可 操作 性 :多 光谱 遥感 影像 的 空间 分 辩 率 
与 高 分 数据 相 比 较 低 ,其 光谱 分 辩 率 与 高 光谱 数据 


的 分 类 体 
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图 1 研究 区 概况 图 


Fig.1 Location of study area 


201903.00158v1 


chinaXiv 


ChinaXiv 合 作 期 刊 


tó wy 


相 比 较 低 ,因此 使 用 多 光谱 遥感 影像 难以 区 分 特征 
相近 的 植被 类 型 。 例 如 , 针 叶 林 和 阔 叶 林 在 多 光谱 
遥感 影像 上 难以 区 分 ,因此 将 他 们 合并 为 乔木 林 。 

参考 《中 国 植 被 》, 将 植被 型 组 作为 分 类 系统 中 
的 最 高 级 别 , 将 具有 相似 生活 型 和 相似 形态 和 外 观 
的 植物 群落 划分 为 相同 的 植被 型 组 。 植 被 型 为 最 重 
要 的 分 类 级 别 ,是 在 植被 型 组 的 基础 上 ,根据 植被 的 
生活 型 和 对 水 热 条 件 的 生态 关系 进一步 划分 的 。 
最 终 建立 的 分 类 体系 如 表 1 所 示 。 


R1 西北 农 牧 交 错 带 植被 分 类 体系 
Tab.1 Vegetation classification system in 


agro-pasturage ecotone of northwest China 


植被 型 组 植被 型 

乔木 林 (A) 温 性 常 绿 针 叶 林 (Al ) 
典型 落叶 阔 叶 林 ( A2) 

灌 从 (C) 常 绿 针 叶 灌 从 (C1) 
温 性 落叶 阔 叶 灌 从 (C2) 

草原 (D) 草 多 草原 (DI1) 
典型 草原 (D2) 
荒漠 草原 (D3) 

荒漠 (下) 荒漠 (El) 

草 向 (H) fa] (H1) 

沼泽 (DT) 沼泽 (1 ) 

农作物 (K) 农作物 (K1) 


1.3 数据 源 

本 文 在 多 光谱 遥感 影像 Landsat8 OLI 数据 的 基 
础 上 ,使 用 地 形 数据 DEM 提供 分 类 所 用 的 辅助 信 
息 。 实 地 考察 和 高 清 有 影像 数据 用 于 辅助 采集 分 类 用 
的 样本 点 。 此 外 ,道路 ,行政 区 划 数 据 等 研究 区 域 的 
基础 数据 主要 用 于 制图 。 

Landsat8 OLI 数 据 下 载 于 USGS 官网 (http: // 
glovis. usgs. gov/) 。 为 了 提取 乔木 林 和 灌 丛 的 常 绿 
和 落叶 信息 ,本 文选 取 研 究 区 内 植被 的 生长 季 和 非 
生长 季 内 质量 较 高 的 影像 各 一 期 ,影像 信息 见 表 2。 
本 文 基于 ENVIS. 1 软件 对 卫星 影像 数据 进行 辐射 
定 标 、 大 气 校正 等 预 处 理 以 获取 真实 地 物 的 反射 率 
信息 。 

30m 的 DEM 数据 下 载 于 地 理 空 间 数据 云 平台 
(http: //www. gscloud. cn/) 。 对 数据 进行 拼接 .投影 
和 裁剪 后 获得 研究 区 的 DEM 数据 。 通 过 对 研究 区 
的 多 次 野外 实地 调查 ,在 积累 了 先 验 知识 的 同时 获 
取 了 大 量 的 样本 点 。 这 些 样本 点 结合 google 高 清 影 
像 可 以 提供 足够 的 样本 点 用 于 影像 分 类 和 精度 验 
证 。 


R2 遥感 影像 信息 表 


Tab.2 Table of image information 


轨道 号 成 像 日 期 云 量 / % 平均 海拔 / m 

127/33 2016/03/13 0.09 1 234 
2015/07/01 0.02 

127/34 2016/02/26 0.09 1172 
2015/07/01 0.01 

128/33 2015/02/14 0.37 1 338 
2016/08/27 0.84 

128/34 2015/02/14 0.27 1 412 
2016/08/27 0.01 

129/32 2015/03/09 0.35 1251 
2015/09/01 0.1 

129/33 2015/03/09 0.25 1 343 
2015/09/01 0.08 

129/34 2016/03/11 0.47 1 483 
2015/09/01 0.27 


14 基于 遥感 的 植被 信息 提取 方法 

西北 农 牧 交错 带 土地 覆 被 情况 较为 复杂 ,使 用 
一 般 的 分 类 方法 难以 取得 较 好 的 分 类 结果 。 因 此 ， 
本 文采 用 ”分 层 分 类 . 逐 层 验证 ”的 方法 提取 植被 信 
息 。 首 先 , 由 于 混合 像 元 分 解 模 型 可 以 有 效 且 准确 
的 提取 植被 覆盖 度 '" , 据 此 可 以 将 遥感 影像 分 为 植 
被 区 域 与 非 植 被 区 域 ;其 次 ,利用 CART 决策 树 结构 
清晰 且 准 确 性 高 ”的 优势 ,将 植被 区 域 进 一 步 分 类 


野外 考察 数据 
数据 预 处 理 ( 师 射 定 标 、 大 气 校正 、 定义 投影 .坐标 转换 、 拼接 、 裁 前 等 


Landsat 遥 感 影像 


相关 指标 获取 
最 佳 指标 筛选 


植被 与 非 植被 分 类 


Ries) 


到 2 技术 路 线 
Fig.2 Technical route of the study 
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为 植被 型 组 ;最 后 ,通过 创建 NDVI 差 值 比值 (NDVI_ 
PR) 指 数 对 类 别 比 较 复杂 的 乔木 林 和 淮 丛 进一步 细 
分 ,得 到 其 落叶 和 常 绿 特征 ,并 通过 使 用 TYD1 值 对 
草原 进行 了 分 类 ,得 到 了 三 种 草原 类 型 ,最 终 获 得 本 
研究 区 的 植被 信息 提取 结果 。 图 2 所 示 为 本 文采 用 
的 技术 路 线 。 


2 植被 信息 的 提取 


2.1 植被 覆盖 区 域 的 提取 

植被 覆盖 度 是 评判 地 表 植 被 覆盖 状况 的 重要 指 
标 之 一 ,因此 也 和 常 被 用 于 区 分 植被 与 非 植被 区 。 本 
文选 取 线 性 光谱 混合 模型 进行 混合 像 元 分 解 以 获取 
研究 区 植被 覆盖 度 信 息 。 该 模型 一 般 包 括 以 下 几 个 
WR: MNF 变换 、 端 元 选取 与 提纯 和 混合 像 元 分 


这 里 使 用 完全 约束 最 小 二 乘法 模型 进行 混合 像 
元 分 解 , 这 一 过 程 可 以 通过 加 载 ENVI 附加 模块 完 
成 。 对 所 得 误差 波段 的 值 进行 统计 ,可 得 大 部 分 像 
元 的 误差 都 较 小 ,以 128/33 为 例 , 其 平均 误差 为 
0.019, 

由 于 研究 区 在 西北 和 干旱 与 半 干 旱 区 ,植被 盖 度 
普遍 偏 低 ,一 般 盖 度 为 5% ~ 20% 。 因 此 ,植被 覆盖 
ERF 5% 的 区 域 均 可 认为 是 有 植被 覆盖 区 域 ( 图 
3) 。 研 究 区 总 面积 为 86 453.44 km ,植被 覆盖 区 域 
总 面积 为 71 450.72 km 。 在 ArcGIS 软件 中 生成 
1 000 个 随机 点 ,并 结合 高 清 影像 判定 提取 结果 的 精 


BE .精度 验证 的 结果 显示 ,提取 结果 总 精度 为 
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Fig.3 Vegetation coverage in the study area 
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93.5% ,Kappa 系数 为 0.782 ,结果 较为 理想 。 
2.2 植被 型 组 分 类 
2.2.1 分 类 指标 ERIRE ,使 用 多 种 
间 标 能 够 有 效 提高 分 类 精度 ,但 是 指标 过 多 也 可 能 
导致 数据 元 余 和 过 分 类 (过 拟 合 ) 。 因 此 ,有 必要 选 
择 合适 的 分 类 指标 进行 分 类 。 不 同 的 地 物 在 各 原始 
波段 有 不 同 的 反射 率 差异 ,但 是 任 一 波段 所 包含 的 
言 息 量 不 同 ,选取 信息 量 最 大 的 几 个 波段 能 够 在 保 
证 分 类 精度 的 同时 提高 分 类 效率 。 在 原始 波段 的 基 
础 上 ,结合 一 些 辅助 信息 可 以 有 效 地 提高 分 类 精度 。 
已 有 的 研究 主要 用 光谱 特征 、 纹 理 特 征 和 地 形 特 征 
完成 植被 的 分 类 。 光 谱 特 征 主要 考虑 植被 冠 层 的 水 
分 含量 .地表 覆 被 的 水 分 含量 以 及 常用 的 多 种 遥感 
分 类 和 环境 监测 指数 ”i 。 本 文选 取 的 光谱 特征 有 
波段 优选 后 的 原始 波段 改进 的 归 一 化 差异 水 体 指 
数 ( Modified Normalized Difference Water Index, 
MNDWI) 归 一 化 差异 湿度 指数 (Normalized Differ- 
ence Moisture Index , NDMI) 以 及 继 帽 变换 、 主 成 分 变 
换 后 的 结果 以 及 植被 指数 NDVI ( Normalized Differ- 
ence Vegetation Index) 。 地 形 特 征 有 数字 高 程 模型 
( Digital Elevation Model, DEM) , 而 纹理 特征 则 选用 
均值 ( 表 3) 。 
研究 表明 '” ,影像 的 原始 波段 的 标准 差 越 大 ， 
相关 系数 越 小 , 则 影像 的 信息 量 越 大 。 以 影像 128/ 
33 为 例 :NIR,SWIR1,SWIR2( 分 别 为 影像 的 第 5 ,6， 
7 波段 ) 的 标准 差 相 对 较 大 ,其 中 ,NIR 与 其 他 各 波 
段 的 相关 系数 均 较 小 ,其 次 为 SWIRI 和 SWIR2 ,而 
两 个 短波 红外 波段 之 间 的 相关 系数 较 大 。 通 过 统计 
7 大 植被 型 组 的 各 600 个 样本 点 ,可 知 不 同 的 植被 
型 组 在 NIR 和 SWIRL 的 差异 相对 较 大 。 因 此 ,本文 
选取 NIR 和 SWIRL 作为 用 于 分 类 的 原始 波段 。 
NDVI?) EVI( 增 强 型 植被 指数 ,Enhanced Veg- 

etation Index) ‘4! SAV ( 土壤 调节 植被 指数 3]， 
Soil-Adjusted Vegetation Index ) , GEMI ( 全球 环 境 监 
测 指 数 , Global Environment Monitoring Index ) [eel 
GRNDVI( 绿 光 与 红 光波 段 归 一 化 植被 指数 ,Green- 
Red NDVI) '*"! 和 TDVI( 转换 型 差 值 植被 指数 ,Trans- 
formed difference vegetation index ) [28] FM H ae af JA 
的 遥感 植被 指数 中 选取 的 适合 干旱 与 半 干 旱 区 的 植 
被 指数 。 通 过 比较 分 析 , 从 这 六 个 指数 中 选取 合适 
的 植被 指数 用 于 分 类 。 定 边 县 北部 的 一 块 区 域 ( 影 
像 128/34 ) 地 表 植 被 覆盖 复杂 ,因此 分 别 在 该 区 域 
计算 上 述 六 个 指标 。 结 果 显 示 NDVI 在 高 植被 覆盖 
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表 3 分 类 指标 


Tab.3 Used indices in the classification 


特征 指标 类 型 分 类 指标 描述 

光谱 特征 原始 波段 NIR ,SWIR1 影像 信息 量 最 大 且 植 被 类 型 光谱 差异 较 大 
植被 指数 NDVI 常用 的 植被 指数 中 分 类 效果 最 好 
水 体 指数 MNDWI 探测 地 表 覆 被 水 分 含量 
湿度 指数 NDMI 研究 植被 冠 层 含水 量 
缕 帽 变换 结果 Brightness , Greenness , Wetness 常用 于 遥感 地 物 分 类 与 遥感 生态 环境 监测 
主 成 分 PC1 ,PC2 ,PC3 减少 影像 元 余 信 息 

纹理 特征 纹理 特征 Mean 各 类 地 物 的 区 分 度 最 大 

地 形 特征 数字 高 程 模 型 DEM 不 同 地 形 条 件 生长 不 同 植被 类 型 


区 和 低 植被 覆盖 区 的 分 类 效果 均 好 于 其 他 几 个。 
此 本 文 最 终 采 用 NDVI 作为 分 类 的 指标 。 

将 纹理 特征 与 光谱 特征 相 结 合 可 以 有 效 提高 遥 
感 分 类 的 精度 ” 。 在 研究 区 内 选取 试验 区 域 (影像 
条 带 号 为 128/34) ,分 别 计算 区 域内 的 8 种 纹理 指 
标 (均值 ,方差 ,对 比 度 , 协 同性 , 相 异 性 , 炉 值 ,二 阶 
BE ,相关 性 ) 。8 种 指标 中 ,均值 方差 和 对 比 度 对 应 
的 标准 差 较 大 ,但 在 均值 的 单 波 段 影像 中 ,不 同 植被 
型 组 呈现 出 不 同 的 色调 , 即 区 分 度 较 大 。 因 此 ,最 终 
选取 均值 (Mean ) 作为 用 于 影像 分 类 的 纹理 特征 。 
2.2.2 植被 型 组 分 类 ”样本 点 选取 。 每 种 植被 类 
型 均 选 取 一 定量 的 样本 点 ,在 选取 样本 点 时 需要 考 
WE (1) 样本 点 要 具有 典型 代表 性 ,选取 样本 点 时 要 
考虑 每 个 植被 类 型 的 各 主要 植物 类 型 ;(2) 样本 点 
应 当 均 匀 分 布 于 整个 研究 区 ; (3 ) 选择 尽 可 能 多 的 
野外 实地 采样 点 ( 共 693 个 点 )。 

统计 汇总 所 有 样本 点 的 所 有 13 个 指标 和 类 别 


AB ,并 将 汇总 结果 导入 Clementine12. 0 软件 中 , 生 
成 分 类 CART 决策 树 。 在 得 到 的 影像 127/33 的 分 
类 决策 树 ( 图 4) 中 可 以 看 出 上 述 13 个 指标 并 未 全 
部 用 到 。 软 件 可 以 自动 生成 评价 决策 树 的 误差 矩 
阵 ,结果 显示 其 标准 中 误差 为 0.003。 在 ENVI 软件 
中 运行 所 得 到 的 CART 决策 树 ,得 到 如 图 5 所 示 的 
植被 型 组 分 类 结果 。 
2.3 植被 型 的 分 类 
2.3.1 乔木 林 和 灌 人 勾 分 类 A201548 月 起 ,项 
目 组 共 对 本 研究 区 进行 四 次 野外 考察 。 在 基本 覆盖 
整个 研究 区 的 多 次 考察 中 ,未 见 到 和 常 绿 针 叶 林 林地 。 
即使 研究 区 确实 有 少量 常 绿 针 叶 林 林地 存在 ,也 很 
难 从 30 m 分 辩 率 的 Landsat 多 光谱 数据 中 将 其 提取 
出 来 。 

在 单一 时 相 的 Landsat 影像 上 区 分 植被 的 落叶 
和 和 常 绿 特征 是 非常 困难 的 ,但 是 由 于 落叶 植被 在 其 
非 生长 季 会 落叶 ,在 影像 上 的 光谱 信息 变化 较为 明 


图 4 基于 多 指标 参与 的 CART 决策 树 


Fig.4 CART decision tree classification model based on multi-index participation 
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到 5 人 研究 区 植被 型 组 分 类 结果 


Fig.5 Classification result of the vegetation type group 


显 。 因 此 ,可 以 使 用 多 时 相 的 影像 来 区 分 落叶 植被 
和 常 绿 植被 。 本 文 基于 植被 生长 季 ( 主要 是 夏季 ) 
和 植被 非 生长 季 ( 主要 是 冬季 ) 的 影像 ,利用 落叶 植 
被 和 常 绿 植被 在 不 同 季节 的 NDVI 的 变化 的 差异 完 
成 两 种 植被 类 型 的 分 类 。 两 种 类 型 的 植被 在 生长 季 
的 NDVI {E ( NDVI_S) 97 FARA KAW NDVI 值 
(NDVI_W) , Ti W 2e te I ZEAE K E PE AE KR EY 
NDVI 值 的 差异 较 小 ,落叶 植被 的 差异 较 大 。 基 于 以 
上 特征 ,创建 了 NDVI 差 值 比值 植被 (NDVI_DR) 指 
数 。 这 一 指数 的 定义 为 :生长 季 的 NDVI 和 非 生长 
ÆW NDVI 的 差 值 和 非 生长 季 的 NDVI 的 比值 , 即 : 


NDVI_S - NDVI_LW 
NDVI_W 


NDVI_DR = (1) 

NDVI_DR 可 以 很 好 地 反映 出 NDVI 在 不 同 季节 
的 变化 强度 。 比 较 了 样本 点 的 NDVI_DR 和 NDVI_ 
形 之 后 ,可 以 得 到 如 下 的 乔木 林 分 类 规则 : 


NDVI_DR > NDVI_W 典型 落叶 阅 叶 林 
| 温带 常 绿 针 叶 林 


通过 统计 常 绿 和 落叶 灌 从 的 样本 点 在 不 同 季 节 
影像 上 的 NDVI 值 , 可知 他 们 与 乔木 林 的 常 绿 和 落 
叶 相 似 , 因 此 可 以 使 用 相同 的 分 类 规则 完成 常 绿 针 
叶 灌 从 和 落叶 治 从 的 分 类 。 

2.3.2 草原 分 类 根据 《中 国 植 被 ) 对 不 同 草原 类 
型 的 定义 可 知 ,不 同类 型 的 草原 有 着 不 同 的 湿度 和 
植被 盖 度 ,因此 区 分 不 同 的 草原 类 型 需要 同时 考虑 
这 两 个 特征 。 温 度 植被 干旱 指数 (Temperature Veg- 
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etation Dryness Index , TVDI) 基于 NDVI 和 地 表 温 度 
很 好 地 结合 了 湿度 和 植被 盖 度 的 特征 。 因 此 本 文通 
过 反 演 得 到 的 7VDI 来 细 分 3 种 草原 类 型 。 在 分 类 
前 ,对 3 种 草原 类 型 各 采集 160 个 样本 点 ,计算 其 
TVDI 值 并 进行 统计 分 析 。 最 终 使 用 不 同 草原 类 型 
的 TV 平均 值 作 为 分 类 冰 值 ,在 ENVI 中 适当 调整 
净值 后 ,建立 草原 分 类 的 决策 树 ,完成 草原 的 精细 化 
分 类 。 最 终 决策 树 使 用 的 分 类 阔 值 为 0. 4091 和 
0.7349 。 

2.4 分 类 后 处 理 及 精度 评价 

分 类 结果 需要 进行 主 次 要 分 析 以 消除 分 类 结 
中 过 于 破碎 的 斑 块 。 在 ENVI 软件 中 , 主 次 要 分 析 
的 变换 核 大 小 有 多 种 。 不 同 的 变换 核对 分 类 结果 的 
影响 程度 不 同 ,其 中 较 大 的 变换 核对 分 类 结果 的 影 
响 较 大 。 为 了 在 优化 分 类 结果 的 同时 减少 这 种 影 
响 ,在 分 析 中 使 用 3 x3 大 小 的 变换 核 。 由 于 所 获取 
的 影像 的 小 部 分 区 域 有 云 覆 盖 , 因 此 这 些 区 域 中 的 
植被 信息 无 法 被 正确 的 提取 。 利 用 其 它 时 相 的 遥感 
影像 作为 辅助 对 分 类 结果 进行 修正 后 ,得 到 最 终 的 
结果 (图 6) 。 

根据 实地 调查 获得 的 先 验 知识 ,基于 google 高 
清 影 像 选 取 用 于 精度 评价 的 样本 。 将 整理 后 的 样本 
点 和 分 类 结果 加 载 到 ENVI 软件 ,计算 混淆 和 矩阵 以 
评价 分 类 结果 的 精度 。 表 4 即 为 计算 所 得 的 精度 评 
价 结果 ( 表 中 0 代表 其 他 , 即 非 植 被 ) 。 最 终 的 分 类 


me 小 性 常 绿 针 叶 林 mm 草包 草原 HWE 
四 典型 落叶 逆 叶 林 mm 典型 草原 mm 沼泽 
at, UL ee 
Se ue OE 农作物 非 植被 区 

106° ÈE 107°E 108° E 109° E 110°E 111° 
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图 6 西北 农 牧 交错 带 植被 信息 提取 最 终结 果 
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表 4 精度 评价 表 
Tab.4 Accuracy evaluation table 

类 型 Al A2 Cl C2 D1 D2 D3 El H1 Il Kl O 总 计 
Al 405 0 42 0 0 0 0 6 2 19 5 3 482 

A2 51 948 30 11 22 2 0 1 3 5 14 1 1 088 
Cl 1 0 194 6 0 9 5 6 9 0 0 0 230 

C2 1 43 60 609 18 17 14 5 17 3 7 0 794 

D1 43 243 0 26 500 37 1 0 42 38 63 11 1 004 
D2 4 19 63 55 38 1 085 56 27 15 21 92 0 1 475 
D3 0 1 0 9 3 44 847 42 1 1 18 0 966 

El 0 0 3 0 0 8 152 817 4 0 3 79 1 066 
H1 5 17 2 2 29 15 0 0 419 65 36 5 595 

ll 3 3 0 0 9 2 0 0 3 350 12 35 417 

Kl 10 18 0 9 27 7 0 0 27 8 1 276 0 1 382 
0 5 0 1 1 4 0 12 31 5 24 9 856 948 
总 计 528 1 292 395 728 650 1 226 1 087 935 547 534 1535 990 10 447 

总 体 分 类 精度 79.51% 


Kappa 系数 


结果 总 体 精 度 为 79. 51% , Kappa 系数 为 0.773 , 相 
较 于 植被 分 布 破碎 复杂 的 实际 情况 ,这 一 分 类 结果 
较 好 。 


3 结论 
本 文 使 用 Landsat 多 光谱 遥感 影像 数据 , 按 "分 


层 分 类 , 逐 层 验 证 " 的 思路 ,对 西北 农 牧 交 错 区 的 植 
被 信息 进行 了 比较 精细 的 提取 ,得 到 以 下 结论 : 

混合 像 元 分 解 模 型 能 够 很 好 地 区 分 植被 和 非 植 
被 的 分 布 范围 。 通 过 在 MNF 变换 后 的 影像 波段 组 
成 的 二 维特 征 空间 中 手动 提取 端 元 ,使 用 n-D Visu- 
alizer 对 所 选 端 元 进行 提纯 ,最 后 使 用 完全 约束 的 最 
小 二 乘 模型 完成 混合 像 元 分 解 ,区 分 有 植被 覆盖 区 
域 和 无 植被 覆盖 区 域 。 验 证 结果 显示 ,有 植被 覆盖 
区 域 的 提取 结果 的 总 体 精度 为 93.5% ,Kappa 系数 
为 0.782 ,结果 较为 理想 。 

基于 多 指标 参与 的 CART 决策 树 能 够 准确 快速 
地 对 复杂 地 表 的 植被 型 组 进行 分 类 。 通 过 筛选 ,使 
用 表征 光谱 特征 、 纹 理 特征 和 地 形 特征 的 13 个 指 
标 , 对 所 选 样 本 点 进行 分 类 ,以 训练 CART 决策 树 。 
最 终 用 于 分 类 的 指标 仅 包 含 NDVI、PC2 SWIR1 
NIR DEM 和 MNDWI 六 个 。 软 件 生 成 的 用 于 评价 
决策 树 的 误差 矩阵 表 显 示 , 其 标准 中 误差 仅 为 0. 
003 ,说 明 决 策 树 的 分 类 结果 较 好 。 

创建 的 NDVI_DR 结合 谱 间 关系 法 能 够 很 好 地 


区 分 植被 的 常 绿 和 落叶 特征 ,使 用 VDI [ed {Et HE 
很 好 地 区 分 不 同 草原 类 型 。 通 过 对 不 同 季 节 NDVI 
的 特征 进行 分 析 ,构建 了 NDVI_DR。 基 于 谱 间 关系 
法 ,构建 了 用 于 区 分 植被 常 绿 和 落叶 特征 的 分 类 规 
则 。 利 用 和 常 绿 植被 和 落叶 植被 在 NDVI_DR 与 NDVI 
_ 丈 大 小 关系 的 差异 ,实现 了 乔木 林 和 灌 丛 的 植被 
型 分 类 。 根 据 草 血 草原 .典型 草原 和 荒漠 草原 在 遥 
感 影 像 上 的 差异 ,使 用 7TVDI 并 设置 阀 值 实现 了 三 
者 的 区 分 。 最 终 分 类 结果 的 总 体 精 度 达 到 79. 
51% ,kappa ABA 0.773. 

本 文 提 出 的 提取 方法 比较 适 于 中 等 分 辨 率 多 光 
谱 双 感 影像 数据 条 件 下 的 大 区 域 信息 提取 ,能 在 提 
高 提取 效率 的 同时 ,减少 主观 因素 的 影响 。 在 以 后 
的 研究 中 ,要 尽 可 能 使 用 多 源 ,多 时 相 的 数据 以 进 一 
步 提高 分 类 精度 。 在 评价 分 类 精度 时 ,应 进一步 降 
低 主观 性 对 评价 结果 的 影响 。 
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Vegetation information extraction in farming-pastoral ectones in northwest 
China using hierarchical classification and multiple indices 
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(1 The college of Earth and Environment , Lanzhou University , Lanzhou 730000 ,Gansu,China; 2 Key Laboratory of Western 
China’s Environmental System (Ministry of Education) ,Lanzhou University , Lanzhou 730000 , Gansu , China) 


Abstract: Based on Chinese vegetation classification criteria described in the book “Chinese Vegetation” , the 
vegetation classification system suitable for the farming-pastoral ectones in Northwestern China was established with 
the field investigation. Using the Landsat images ,the terrain data and field data of the study area, the vegetation in- 
formation in the study area was refined according to the implementation strategy as “hierarchical classification and 
layer by layer verification”. During the extraction process ,the mixed pixels decomposition of preprocessed remotely 
sensed images was performed by using the fully constrained least-squares model, and the vegetation fraction of the 
study area was obtained. We classified the study area into vegetation area where the vegetation fraction is larger than 
5% and non-vegetation area where the vegetation fraction is less than 5% . In the vegetation area ,it was further clas- 
sified into 7 main vegetation type groups which include the tree group,shrub group and grassland group using the 
CART decision tree based on the spectral characteristic , texture characteristic (Mean) and terrain characteristic 
( Digital Elevation Model, DEM). Each vegetation type group was again further classified into different sub types 
based on the refined indexes. The tree group and shrub group were categorized into evergreen vegetation type and 
deciduous vegetation type based on the NDVI difference ratio index which was established using the seasonal varia- 
tions of their NDVIs of different plants. The grassland type group was categorized into desert grassland, the typical 
grassland and meadow grassland using the temperature vegetation dryness index ( TVDI/). After this step ,the spatial 
distribution of each vegetation type was obtained. It was proved that the overall accuracy of the final classification 
can reach 79.51% and the kappa coefficient is 0.773. The classification method used in this study makes full use 
of the spectral information and texture information of the remotely sensed images ,and cooperates with terrain infor- 
mation. Experimental result shows that the method of using hierarchical classification and multiple indexes could ex- 
tract the vegetation information efficiently from the images of the farming-pastoral ectones with high accuracy. The 
classification result in this area provides the basic data for the further research on the relationship between surface 
hydrothermal process and land cover change , especially vegetation cover change. Meanwhile ,it provides reference for 
the conservation of vegetation area and ecological environment construction in this area. 

Key words: the farming-pastoral ectone; vegetation information extraction; LSMA; CART decision tree; NDVI 

difference ratio (NDVI_DR) ; temperature vegetation dryness index ( TVDI) 


